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ABSTRACT 

In this work, considering the impact of a supernova remnant (SNR) with a neutral 
magnetized cloud we derived analytically a set of conditions that are favorable for 
driving gravitational instability in the cloud and thus star formation. Using these 
conditions, we have built diagrams of the SNR radius, Rsnb., versus the initial cloud 
density, ric, that constrain a domain in the parameter space where star formation is 
allowed. This work is an extension to previous study performed without considering 
magnetic fields (Melioli et al. 2006). The diagrams are also tested with fully 3-D MHD 
radiative cooling simulations involving a SNR and a self-gravitating cloud and we find 
that the numerical analysis is consistent with the results predicted by the diagrams. 
While the inclusion of a homogeneous magnetic field approximately perpendicular to 
the impact velocity of the SNR with an intensity ~ 1 fiG within the cloud results 
only a small shrinking of the star formation zone in the diagram relative to that 
without magnetic field, a larger magnetic field (~ 10 fiG) causes a significant shrinking, 
as expected. Though derived from simple analytical considerations these diagrams 
provide a useful tool for identifying sites where star formation could be triggered by 
the impact of a SN blast wave. Applications of them to a few regions of our own galaxy 
(e.g., the large CO shell in the direction of Cassiopeia, and the Edge Cloud 2 in the 
direction of the Scorpions constellation) have revealed that star formation in those sites 
could have been triggered by shock waves from SNRs for specific values of the initial 
neutral cloud density and the SNR radius. Finally, we have evaluated the effective star 
formation efficiency for this sort of interaction and found that it is generally smaller 
than the observed values in our own Galaxy (sfe ~ 0.01—0.3). This result is consistent 
with previous work in the literature and also suggests that the mechanism presently 
investigated, though very powerful to drive structure formation, supersonic turbulence 
and eventually, local star formation, does not seem to be sufficient to drive global star 
formation in normal star forming galaxies, not even when the magnetic field in the 
neutral clouds is neglected. 
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Essentially all present day star formation takes place in 
molecular clouds (MCs; e.g. Blitz 1993; Williams, Blitz & 
McKee 2000). It is likely that the MCs are relatively tran- 
sient, dynamically evolving structures produced by compres- 
sive motions in the diffuse HI medium of either gravitational 
or turbulent origin, or some combination of both (e.g., Hart- 
mann et al. 2001; Ballestero-Paredes et al. 2007). In fact. 
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observations of supersonic line-widths in the MCs support 
the presence of supersonic turbulence in these clouds with a 
wealth of structures on all length-scales (Larson 1981; Blitz 
& Williams 1999; Elmegreen & Scalo 2004; Lazarian & Es- 
quivel 2003). Recent numerical simulations in periodic boxes 
have shed some light on the role of the turbulence in the evo- 
lution of the MCs and the star formation within them. They 
suggest that the continuous injection of supersonic motions, 
maintained by internal or external driving mechanisms (see 
below), can support a cloud globally against gravitational 
collapse so that the net effect of turbulence seems to inhibit 
collapse and this would explain the observed low overall star 
formation efficiencies in the Galaxy (Klessen et al. 2000; Mac 
Low & Klessen 2004; Vazquez- Semandeni et al. 2005). On 
the other hand, the supersonic turbulence is also able to 
produce density enhancements in the gas that may allow 
local collapse into stars in both nonmagnetized (Klessen et 
al. 2000; Elmegreen & Scalo 2004) and magnetized media 
(Heitsch et al. 2001; Nakamura & Li 2005). 

What are the possible sources of this turbulence in the 
MCs? Suggested candidates for an internal driving mecha- 
nism of turbulence include feedback from both low-mass and 
massive stars. These later, in particular, are major structur- 
ing agents in the ISM in general (McCray & Snow 1979), 
initially through the production of powerful winds and in- 
tense ionizing radiation and, at the end of their lives through 
the explosions as supernovae (SNe). It is worth noting how- 
ever, that MCs with and without star formation have sim- 
ilar kinematic properties (Williams, Blitz & McKee 2000). 
External candidates include galactic spiral shocks (Roberts 
1969; Bonnell et al. 2006) and again SNe shocks (Wada & 
Norman 2001; Elmegreen & Scalo 2004). These processes 
seem to have sufficient energy to explain the kinematics of 
the ISM and can generate the observed velocity dispersion- 
sizescale relation (Kornreich & Scalo 2000). Other mecha- 
nisms, such as magnetorotational instabilities, and even the 
expansion of HII regions and fluctuations in the ultraviolet 
(UV) field apparently inject energy into the ambient medium 
at a rate which is about an order of magnitude lower than 
the energy that is required to explain the random motions 
of the ISM at several scales, nonetheless the relative impor- 
tance of all these injection mechanisms upon star formation 
is still a matter of debate (see, e.g., Joung & Mac Low 2006; 
Ballesteros-Paredes et al. 2007; Mac Low 2008 for reviews). 

In this work, we focus on one of these driving mech- 
anisms — supernova explosions produce large blast waves 
that strike the interstellar clouds, compressing and some- 
times destroying them, but the compression by the shock 
may also trigger local star formation (Elmegreen & Lada 
1977; see also Nakamura & Li 2005 and references therein). 
There has been extensive analytical and numerical work ex- 
ploring the role of the SNe in generating supersonic tur- 
bulence and multi-phase structures in the ISM (e.g., Cox 
& Smith 1974; McKee & Cowie 1975; McKee & Ostriker 
1977; Kornreich & Scalo 2000; Mac Low & Klessen 2004; 
Vazquez-Semadeni et al. 2000; de Avillez & Breitschwerdt 
2005) and even at the galactic scale, in the generation of 
large-scale outflows like galactic fountains and winds (e.g., 
de Gouveia Dal Pino & Medina-Tanco 1999; de Avillez 2000; 
Heckman et al. 2001; de Avillez & Berry 2001; Melioli et al. 
2008a, 2008b). The coUective effect of the SNe is hkely to be 
the dominant contributor to the observed supersonic turbu- 



lence in the Galaxy (Norman & Ferrara 1996; Mac Low & 
Klessen 2004), however, it is possible that it inhibits global 
star formation rather than triggering it (Joung & Mac Low 
2006). Here, instead of examining the global effects of mul- 
tiple SNe explosions upon the evolution of the ISM and the 
MCs, we will explore the local effects of these interactions. 
To this aim, we will consider a supernova remnant (SNR) 
either in its adiabatic or in its radiative phase impacting 
with an initially homogeneous diffuse neutral cloud and show 
that it is possible to derive analytically a set of conditions 
that can constrain a domain in the relevant parameter space 
where these interactions may lead to the formation of grav- 
itationally unstable, collapsing structures. In Melioli et al. 
(2006; hereafter Paper I), we gave a first step into this di- 
rection by examining interactions involving a SNR and a 
non-magnetized cloud and then applying the results to a SF 
region in the local ISM where the young stellar association of 
P Pictoris was born. Presently, besides including the effects 
of the magnetic fields, we will extend this analysis to few 
other SF regions with some indication of recent-past inter- 
actions with SN shock fronts (e.g., the Edge Cloud 2, Yasui 
et al. 2006, Ruffle et al. 2007; and the "Great CO Shell", 
Reynoso & Mangum 2001). Finally, we will also test our an- 
alytically derived SF domain with 3D MHD simulations of 
SNR interactions with self-gravitating neutral clouds. 

In Section [2l we consider the equations that describe 
the interaction between an expanding SNR and a cloud. In 
Section [31 we obtain analytically the set of constraints from 
these interactions with both non-magnetized and magne- 
tized clouds that may lead to star formation (SF) and then 
build diagrams where these constraints delineate a domain 
in the parameter space which is appropriate for SF. In Sec- 
tion|4l we describe radiative cooling 3D MHD simulations of 
these interactions which confirm the results obtained from 
the analytical study, and[5]we apply these analytical results 
to few examples of regions of the ISM that present some 
indication of SF due to a recent SNR-cloud interaction. In 
Section [6l we briefly discuss the role of these interactions to 
the global efficiency of star formation in the Galaxy and in 
Section [3 we draw our conclusions. 



2 SNR-CLOUD INTERACTIONS: AN 
ANALYTICAL DESCRIPTION 

As in paper I, we start by considering the equations that are 
relevant in the study of the interaction between a SNR and 
a cloud. A type II SN explosion generates a spherical shock 
wave that sweeps the interstellar medium (ISM), leading to 
the formation of a SNR. The interaction between a SNR 
and a cloud may compress the gas sufficiently to drive the 
collapse of the cloud. In order to describe analytically this 
interaction we will consider a diffuse neutral cloud with ini- 
tially homogeneous density and constant temperature. After 
the impact, an internal forward shock propagates into the 
cloud with a velocity Vcs- The ram pressure of the blast 
wave, ~ nshVgf^fj, must be comparable to the ram pressure 
behind the shock in the cloud, ~ ncuL and this results the 
well known relation for Vcs in the case of a planar shock: 

v..~vsNn[—) (1) 
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where (^^) = X is the density contrast between the SNR 
shell and the cloud, riah is the shell density and tie is the 
gas cloud density. During its evolution a SNR undergoes 
two main regimes: an adiabatic (or Sedov- Taylor) and later 
a radiative regime. When the SNR is still in the adiabatic 
phase, we may use its expansion velocity (as derived in Eq. 
4 of Paper I) to find 
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where ric, lo is the cloud density in units of 10 cm~ , E^i 
is the SN energy in units of 10^^ erg, Rsnr,50 is the SNR 
shell radius in units of 50 pc. When the SNR enters the 
the radiative phase, the expansion velocity of the SNR is 
described by Eq. (7) of Paper I and this results a forward 
shock into the cloud: 
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where /lo (/ = (^^)) is the density contrast between the 
SNR shell and the ISM density in units of 10 and n is the 
ambient medium density. 

The equations above are valid for a planar shock only. 
Since it is not rare that both the SNR and the cloud have 
dimensions which are of the same order, the effects of cur- 
vature of the shock should be considered. In Ap.|Al we have 
derived a correction factor (I) that must be multiplied to Ucs 
for the case of spherical cloud— SNR interactions, i.e.: 

(4) 

where I is given by eq. HAlf) of Ap. |^ 



2.1 The inclusion of the magnetic field in the 
cloud 

Zeeman measurements indicate that the magnetic flux-to- 
mass ratios in molecular cloud cores are close to the critical 
value for magnetic support (e.g., Crutcher 1999, 2005, 2008), 
while in neutral clouds these ratios are observed to be su- 
percritical. In other words, the typical observed magnetic 
fields in these clouds (~ 6fiG - lO/iG) are smaller than the 
critical value {Bcr — 2nG^^^'N, where N is the cloud sur- 
face density) necessary to suppress gravitational instability 
in them (e.g., Nakano & Nakamura 1978). This could be an 
indication that the presence of the magnetic field in these 
neutral clouds is not relevant when considering their inter- 
actions with SNRs. However, the impact will compress the 
magnetic field behind the shock and this may affect the evo- 
lution of the shocked material and therefore, its conditions 
for gravitational collapse. 

Let us then consider a magnetized diffuse cloud. For 
simplicity, let us assume that its magnetic field is initially 
approximately uniform and, in order to maximize the effects 
of the field we shall also assume that it is initially normal to 



the SNR shock velocity at the impact. The observed mag- 
netic fields are actually randomly distributed. This means 
that only a fraction of its average strength will be effectively 
normal to the SNR velocity and work against the impact. 
We will find below that when we take an effective value for 
the cloud magnetic field Be ~ 1/iG, its effect upon the SF 
diagrams is in fact not too relevant. However, when we take 
values ~ 5 — 10 ^G, the allowed SF domain in the diagrams 
shrank considerably. 

Depending on the physical conditions of the 
cloud, the propagation of the shock front into it can 
be either adiabatic or radiative. For the shocked gas 
at temperatures T ^ lO'' K, we find that the radia- 
tive cooling time is shorter than the crushing time 
(see e.g. Melioli, de Gouveia Dal Pino &: Raga 2005) 
and therefore, we can assume that the forward shock 
wave propagating into the cloud is approximately 
radiative. The Rankine-Hugoniot relations for a ra- 
diative strong shock (with M > 10) (see, e.g. Draine 
& McKee 1993), are 



Tc^sh — Tc 
nc,sh,B ^ y Tie 

Bc,sh ~ y Be 



y = 



2M-2 + MX + [(2M-2 + M^y + 8M-^]i/2 



(5) 
(6) 
(7) 

(8) 



where Tc^sh, nc^sh,B and -Bc.sh are the temperature, density, 
and magnetic field of the shocked cloud gas, respectively, 
M is the Mach number given by Eqs. (|A8|I and ()A9|I for an 
adiabatic and for a radiative-phase SNR, respectively, and 
Ma is the Afvenic Mach number. For an interaction with a 
SNR in the adiabatic phase, using equation ([2]), Ma becomes 
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where va is the Alfven speed in the cloud, _Bc,6 is the pre- 
shock magnetic field of the cloud in units of 10~^ G and 
is the factor I (eg. lAip calculated for Rsnr/tc ~ 5. For an 
interaction with a SNR in the radiative phase, using Eq. 
Ma is given by: 



MA.r = 75 
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(10) 



Replacing (|9]) and (|A8|) in Eq. (|6]), we obtain the density of 
the shocked gas in the magnetized cloud after the interaction 
with an adiabatic SNR: 



nc,sh,B,a 



where 
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Using (|10|) and (|A9|) . we obtain the density of the shocked 
cloud gas after the interaction with a radiative SNR: 
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where 

H — 1 Tc, 100 n-c,io + 0.18 -Bc,6 
and 
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When 1 ^ Ma ^ Af^, the equations above become 
(e.g., Draine & McKee 1993): 



nc,sh,a = 2^''^Ma ric ~ 970 ric 



and 
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(13) 



(14) 



In the limit 1 <^ <^ Ma, which imphes a dynamically 
negligible magnetic field behind the shock, we recover the 
equations without magnetic field derived in Ap. 1X1 fEg. lAlOl 
and lAllI see also Paper I). 



3 CONDITIONS TO DRIVE GRAVITATIONAL 
COLLAPSE IN THE CLOUD 

We can now derive a set of constraints that will determine 
the conditions for driving gravitational instability in the 
cloud right after the interaction with a SNR. A first con- 
straint will simply determine the Jeans mass limit for the 
compressed cloud material. A second one will establish the 
condition on the SNR forward shock front at which it will 
not be too strong to destroy the cloud completely making 
the gas to disperse in the interstellar medium before becom- 
ing gravitationally unstable. A third constraint will estab- 
lish the penetration extent of the same shock front inside 
the cloud before being stalled due to radiative losses. The 
shock must have energy enough to compress as much cloud 
material as possible before fading. Let us derive these three 
constraints. 



3.1 The Jeans Mass Constraint 

3.1.1 In the absence of magnetic field 

In Paper I, we have derived the Jeans mass limit for the 
shocked gas in interactions involving non-magnetized clouds. 
The more precise equations obtained in Ap.|A]for spherical 
SNR— cloud interactions produce slight modifications in the 
Jeans mass derived in Eqs. (22) and (23) of paper I. When 
the interacting SNR is still in the adiabatic regime, the cor- 
responding Jeans mass of the shocked cloud material in the 
absence of magnetic field is given by: 
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and if the SNR is in the radiative regime: 
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In terms of the SNR radius, the conditions above 
over the shocked gas of the cloud imply: 

j2/3 2/3 2 

HsNR.a <, 55.8 ^j^g pc 



t: 



and 

RsNR,r < 55.1 

respectively. 
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(18) 



3.1.2 In the presence of magnetic field 

When including the magnetic field in the cloud, the cor- 
responding minimum mass (Jeans mass) that the shocked 
material must have in order to suffer gravitational collapse 
is given using the Virial Theorem by 
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In terms of the cloud pre-shock and the SNR parameters 
this condition reads (using Eqs. [5]- [Tjl: 
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We can obtain Rsnr as function of Uc and then obtain 
an approximate condition for gravitational collapse solving 
the condition lUc ^ mj,B, where nic is the cloud mass, or: 

,1/3 



: 10 y > 5.63 ric, 10 T'c,ioo + 0.54 B^e y 



(21) 



where y is given by Eq. ([8]). Replacing y in the equation 
above and substituting Eqs. HA8|) and (|9]) we obtain numer- 
ically the Jeans limit or an upper limit for RsNR,a for the 
shocked gas of a magnetized cloud due to the impact with 
a SNR in the adiabatic regime and substituting (|A9|I and 
(|10|) we obtain the same conditions for an interaction with 
a SNR in the radiative regime (see Section 3.4 below). 

We notice that in the limit that ^ p the equa- 
tion above (|2H) recovers the solutions (|17p and (|18p . for the 
adiabatic and radiative phases, respectively. 



3.2 Constraint for non-destruction of the cloud 
due to a SNR impact 

3.2.1 In the absence of magnetic field 

As remarked before, if the SNR-cloud interaction is too 
strong this may cause the complete destruction of the cloud 
even if the shocked material had a priori a total mass 
superior to the Jeans limit. The stability of a cloud 
against destruction right after the impact of a wind 
or a SNR due to the growth of Rayleigh- Taylor and 
Kelvin-Helmholtz instabilities has been explored in 
detail by several authors (see e.g., Murray et al. 
1993; Poludnenko et al. 2002; Melioli, de Gouveia 
Dal Pino & Raga 2005, and references therein). In 
order to obtain the condition for a shocked cloud not 
to be destroyed, we must compare the gravitational 
free-fall time-scale of the shocked gas with the de- 
struction time-scale due to the impact and the con- 
sequent development of those instabilities, td. More 
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precisely, in order to have collapse, a gravitationally unsta- 
ble mode (with typical time tun) must grow fast enough 
to become non-linear within the time-scale of the cloud- 
SNR shock interaction (see, e.g., Nakamura et al. 2006). 
Following Nakamura et al., the condition for non- 
destruction then reads: tun < td, where td is a few 
times the crushing time, tec = 2rc/vc3 (eqs. A3 and 
A4), i.e., the time that the internal forward shock 
takes to cross the cloud. With the help of numerical 
simulations that have taken into account the effects 
of non-equilibrium radiative cooling in the interac- 
tions between clouds and SNRs, Melioli, de Gouveia 
Dal Pino &: Raga (2005) have shown that the cloud 
destruction time is td ~ 4-6 tec, when radiative cool- 
ing is present. Hence, in order to have collapse, we should 
expect that tun < 5tcc- This implies a Mach number: 



where 



M < 60.5 



7/9 



(22) 



where ~ 0.8 is the shock compression factor in the 
transverse direction (Nakamura et al. 2006) and 

{G'pcy^' 



fio = rUc 



(23) 



is the non-dimensional cloud mass, and mc is the 
cloud mass. Thus, we have: 
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We note that in Paper I, the exponent of Vc in the equation 
above was incorrectly given by 3. This relation implies: 
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for an interaction with a SNR in the adiabatic regime, and 

TpO.33 rO.2 rpO.26 rO.4 

RsNR,r > 48.7 0.7 „0.17^0.93 (26) 

'V,10 " 'c.lO 

for an interaction with a SNR in the radiative phase. 



3.2.2 In the presence of magnetic field 

When the cloud is magnetized with an average magnetic 
field normal to the shock velocity then the shock strength 
must be decreased by the magnetic pressure and tension 
forces. Using the same condition for collapse as described in 
the previous paragraph, but including the magnetic field in 
the cloud, we obtain the following condition over the Mach 
number: 
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(27) 



1 We notice that in paper I, we have assumed a condition for 
collapse which was given by tun < 3tcc. In fact, this is a more 
appropriate condition for shocks where the radiative cooling is 
not much important. Presently, we have corrected this condition 
by tun < 5tcc, as td ~ 5tcc is a more representative average 
value for the cloud destruction time in the presence of a strong 
radiative cooling shocked cloud (Melioli, de Gouveia Dal Pino & 
Raga 2005). 
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Substituting the relations found for y (Eq. [S}, M and Ma 
for a collision with an adiabatic (Eqs. |9] and IA8|) and with 
a radiative SNR (Eqs. [TUl and IA9] ) into the equation above 
(Eq. I27|l we can find numerically the new constraints over 
RsNR,a and RsNR,r, respectively, in order to not destroy the 
magnetized cloud at the impact and allow its gravitational 
collapse (see below). 



3.3 Penetration extent of the SNR shock front 
into the cloud 

Besides the constraints derived in the previous sections (3.1 
and 3.2), we still need a third condition upon the impact. 
As remarked before, the shock should have strength enough 
to penetrate into the cloud and compress as much material 
as possible for this to undergo gravitational collapse. 



3. 3. 1 In the absence of magnetic field 

When the shock propagates into the cloud, it will 
decelerate due to continuous radiative losses and 
to the resistance of the unshocked cloud material. 
We can estimate the approximate time at which the 
shock will stall within the cloud (tst) using energy 
conservation arguments. At t = tst, the velocity of 
the shocked gas in the cloud should become ~ Cs, 
so that the Mach number will decay from its initial 
value right after the impact to M ~ 1. As in paper I, 
from the energy conservation behind the shock, we 
find approximately that: 



2 ixmn 2 



5 fes Tc,sh{tst) 1 2 \ , ^\I'c,sh{tst)] nc,sh{tat) , 

— o ^ o ""shV-st) H tat {-ia) 

Z t-tmn 2. finriH 

where the left-hand side is the total energy behind 
the shock in the cloud immediately after the impact 
and the right-hand side is the total energy behind 
the shock at the time tat when it stalls. The ini- 
tial temperature and density of the shocked cloud 
material in the left hand-side of the equation are 
approximately given by the adiabatic values. On 
the right hand side, at tat, the shocked tempera- 
ture and density are obtained from the Rankine- 
Hugoniot relations for a radiative shock with Af ~ 1, 
i.e., Tc,ah{tat) — Tc, and nc,ah{tat) — nc. The radiative 
cooling function at tat, ^[Tah{tat)] — A{Tc), can be ap- 
proximated by that of an optically thin gas (see e.g. 
Dalgarno & McCray 1972). The substitution of these 
conditions into the equation above results: 



We notice that in Paper I, there was a typo in the equation 
for tat (Eq. (31) of that paper). The multiplying factor in that 
equation (9/16) has been now correctly replaced by 31/32. 
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31 



■Vcs 



32 nck{Tc) 



(29) 



where u^s is given by eq. (|4|) or, in terms of the ini- 
tial Mach number, by eq. (|A8|) or (|A9p . and A(Tc) is 
calculated for Tc = 100 K. 

As in paper I, we use this time to compute the maximum 
distance that the shock front (initiated by a given SNR) can 
travel into the cloud before being stopped. This distance is 
then compared with the diameter of the cloud (vcs tst ^ 
2rc) in order to establish the maximum size (that is, the 
minimum energy) that the SNR should have in order to 
generate a shock able to compress the cloud as much as 
possible before being stalled. In the case of an adiabatic 
SNR this gives: 
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And in the case of a SNR in the radiative regime, we find: 
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where A27 is the cooling function (A) in units of 10~^^ erg 
cm"^ s~^, which is the approximate value for Tc — 100 
K and an ionization fraction 10^^ (see below). 



3.3.2 In the presence of magnetic field 

When considering the presence of the normal magnetic field 
in the cloud, our estimates from energy conservation imply 
approximately : 
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(32) 



Where again the left-hand side gives the total en- 
ergy behind the shock in the cloud immediately after 
the impact with the initial temperature, density and 
magnetic field of the shocked cloud material being 
approximately given by the adiabatic values. The 
right-hand side gives the total energy behind the 
shock at tst when Vsh decreases to Vsh{tstY — ICs+u^- 
At this time, Uc^shitst), T^^shitst), and Bc,sh(tst) are ob- 
tained from the Rankine-Hugoniot relations for a 
magnetized radiative shock (eqs. [5]to[7|), for y ~ 1. 
So that again A[Tshitst)] ^ A(rc). The substitution of 
these conditions into the equation above results approxi- 
mately the same estimate for tst as in Eq. (|29|) for the non- 
magnetic case and, therefore, the same upper limits over 
RsNR in order to the shock to sweep the cloud as much as 
possible before being stalled (Eqs.[30landl31[). This is due to 
the fact that for the typical fields observed in these clouds 
(~ 10"* - lO"** G), the magnetic energy terms in Eq. ((32)) 
are negligible compared to the others. 



3.4 Diagrams for Star Formation 

The three constraints derived above in Sections 3.1, 3.2, and 
3.3 for both non-magnetized and magnetized clouds inter- 
acting with SNRs either in the adiabatic or in the radiative 
phase can be plotted together in a diagram showing the SNR 
radius versus the initial (un-shocked) cloud density for dif- 
ferent values of the cloud radius, as performed in Paper I. 
Figures 1 and 2 show examples of these diagrams for non- 
magnetized and magnetized clouds with B — 1 fiG, respec- 
tively, with an initial temperature Tc = 100 K and radius 
varying between rc = 1 pc and 20 pclf| The three constraints 
establish a shaded zone in the parameter space of the di- 
agrams where conditions are appropriate for gravitational 
collapse of the shocked cloud material. Only cloud— SNR 
interactions with initial physical conditions (rc, ric, and 
RsNR)lyir^g within the shaded region (between the solid, 
dotted and dashed lines) may lead to a process of star for- 
mation. We have chosen a SNR in the adiabatic phase in 
these diagrams because it has more stored energy than one 
with the same characteristics in the radiative phase (see, 
however. Figure]^ an example of an interaction involving a 
radiative SNR). As in paper I, we should remark that ac- 
cording to the equations ^ [TH [25l [26l |30] and |3T}, these 
interactions are not very sensitive to the initial temperature 
of the cloud and this explains why we have taken only a char- 
acteristic value for it. We further notice that a cloud with 
a temperature in the range of 10—50 K and a radius larger 
than 10 pc is already Jeans unstable over a large range of 
densities ( > 5 cm~^ when the magnetic field is neglected) 
and does not require, in principle, a shock wave to trigger 
SF. Besides, it will be more difficulty for a SNR shock front 
to destroy a cloud at these temperatures. 

Figure 1, which describes interactions with a non- 
magnetized cloud, was already presented in Paper I. How- 
ever, the modifications in the equations described in Sections 
13. 11 - 13. 3l above have resulted slight modifications in the dia- 
grams. Few remarks are in order with regard to this Figure: 



(i) In paper I, the dotted (blue) curves of the diagrams 
were built for one value only of the radiative cooling func- 
tion of the shocked material, i.e. A{Tsh) — 10~^^ erg cm"^ s~^ 
which is valid for a diffuse gas with temperature 100 K and 
ionization fraction ^ 10~*. Considering that the constraint 
established by the dotted (blue) curve is sensitive to A{Tsh) 
(through Eq. I30p which in turn, can vary by two orders of 
magnitude depending on the value of the ionization fraction 
of the cloud gas, we have presently plotted in the diagrams 
three different dotted (blue) curves in order to cover a more 
realistic range of possible ionization fractions from 0.1 to 
10~*, corresponding to A{Tsh) = lO"'^^ erg cm'' s~^ (lower 
dotted curve), 5 x 10~^^ erg cm"' s~^ (middle dotted curve), 
and 3 x 10~^^ erg cm'' s~^ (upper dotted curve), respec- 
tively (see Dalgarno & McCray 1972). The middle dotted 
(blue) curve corresponds to the average value of A in the 



We notice that in the diagrams of Figs. 1 and 2, and in 
most of the diagrams of this work, we have assumed an 
ambient density n = 0.05 cm^^. We have chosen this low 
density value in order to try to better reproduce the hot 
phase, low density medium commonly expected around 
a SNR, particularly in its adiabatic phase. 
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Figure 1. Constraints on the SNR radius versus cloud density for 
4 different cloud radius. Prom top to bottom: first panel, rc = 1 pc; 
second panel, rc = 5 pc; third panel, rc = 10 pc; and forth panel, 
rc = 20 pc. Dashed (green) line: upper limit for complete cloud de- 
struction after an encounter with an adiabatic SNR derived from 
Eq. II25I I; solid (red) line: upper limit for the shocked cloud to 
reach the Jeans mass derived from Eq. I|17| l: dotted (blue) lines: 
upper limits for the shock front to travel into the cloud before 
being decelerated to subsonic velocities derived from Eq. I|30|l for 
different values of the cooling function A{T^fi) = lO"'^^ erg cm'^ 
s~^ (lower curve), 5 X 10~^^ erg cm^ (middle curve), and 
3 X 10~^^ erg cm^ s~^ (upper curve). The shaded area (between 
the solid, dashed and dotted lines) defines the region where star 
formation can be induced by a SNR-cloud interaction. The crosses 
and the triangle in the panels indicate the initial conditions as- 
sumed for the clouds in the numerical simulations described in 
Section 4.2 of Paper I (see text for details). 



Local Star formation triggered by SN 7 

range above, 5 x 10^^^ erg cm'^ s~^, and could be taken as 
a reference. 

(ii) In the solution presented in Paper I for the cloud with 
Vc = 1 pc (top panel of Figure [T]), there was no allowed 
SF shaded zone. According to the present corrections and 
modifications, we see that a thin shaded "star-formation un- 
stable" zone appears now when the cooling function A has 
values which are smaller than 10~^'' erg cm'^ s~^, or ioniza- 
tion fractions <; 10~^. 

(iii) The cross labeled in the bottom panel of Figure 1 
for an rc = 20 pc diffuse cloud corresponds to the initial 
conditions of the numerical simulations presented in Figure 
6 of paper I (i.e., for a SNR at a distance Rsnr ~ 42 pc 
from the surface of the cloud). In the paper I, that cross 
lies outside the unstable shaded zone just above the upper 
limit for a complete shock penetration into the cloud (the 
middle dotted, blue line of the diagram). With the present 
modifications, the cross now lies near the upper limit of 
the shaded unstable zone for values of the cooling function 
A < 10^^^ erg cm^ s^^, or ionization fractions < 10~^. This 
result remarks how sensitive the analytical diagrams are to 
the choice of A (or the ionization fraction) for a given ini- 
tial cloud temperature. According to the radiative cooling 
chemo-hydrodynamical simulations of Figure 6 of Paper I 
(which corresponds to the cross in the diagram of Fig. [1} , the 
SNR shock front really stalls within the cloud before being 
able to cross it completely and, after all the shocked cloud 
material does not reach the conditions to become Jeans un- 
stable, as predicted, but the dense cold shell that develops 
may fragment and later generate dense cores, as suggested 
in Paper I. This points to an ambiguity of the results due to 
their sensitivity to A and the real ionization fraction state 
of the gas. We should also remark that the constraint estab- 
lished by the dotted (blue) curves in the diagrams is actually 
only an upper limit for the condition of penetration of the 
shock into the cloud. The condition that the mach number 
of the shocked material goes to M ~ 1 implies the pressure 
equilibrium between the shocked and the unshocked cloud 
material at the time tst when the shock stalls within the 
cloud. A quick exam of the numerical simulations of Fig- 
ure 6 of Paper I, however, shows that the shock front stalls 
even before this balance is attained. This implies that the 
time tst could be possibly smaller and therefore, the dotted 
(blue) curves in the diagrams should lie below the location 
predicted by Eq. (|30} . 

In spite of the important alterations above in the dia- 
grams of Figure [T] the main results and conclusions of Pa- 
per I for interactions involving SNRs with non — magnetized 
clouds, particularly those regarding the young stellar system 
/3-Pictoris, have remained unchanged (see below, however, 
the implications for this system when the magnetic field is 
incorporated into the cloud). 

Figure [5] shows the same diagrams as in Figure [T] ex- 
cept that they include the effects of the magnetic field in 
the cloud, as described by Eqs. 1211 1271 and 1301 We notice 
that the presence of a normal magnetic field to the shock 
front with an intensity of 1 inhibits slightly the do- 
main of SF in the diagrams, as expected. The magnetic field 
plays a dominant role over the Jeans constraint (the solid 
(red) line in the diagrams) that causes a drift of the 
allowed (shaded) zone of SF to the right in the diagrams 
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(i.e., to larger cloud densities) when compared to the dia- 
grams without magnetic fields of Fig. [1] This drift how- 
ever must be interpreted with care. When deriving 
the Jeans constraint in the presence of a non-null 
uniform magnetic field normal to the shock front 
(eq. I19p we assumed an one-dimensional collapse. 
However, although the presence of B may affect the 
initial compression and collapse of the shocked ma- 
terial, the later evolution and collapse of this ma- 
terial in three-dimensional space will occur mostly 
in the direction parallel to B and thus will not be 
any further affected by it. For this reason, we should 
expect that the actual domain of an unstable magne- 
tized cloud in the diagrams of Figure 2 will be more 
extended to the left of the solid (red) line and will 
be ultimately bounded by the dotted-dashed (red) 
line that gives the Jeans constraint for a null (or 
parallel) magnetic field, just like in Figure 1. This 
part of the diagrams with non-null B will be partic- 
ularly important when comparing them with the 3- 
D numerical simulations of cloud-SNR interactions. 
When larger magnetic field intensities are considered (5-10 
^G) there is a significant shrinking of the allowed SF zone 
in the diagrams (see Figure [3} . This can be understood in 
terms of the mass to magnetic fiux ratio of the cloud before 
the impact. It is given by N/B — 885(nc,io''c,io)/Bc,6. This 
implies B/Bcr ~ Q-TBc^ / {nc,\o'rc,io), which is larger than 
1 for nc,io 5; 34.8 cm~'^ and nc,io ^ 69.6 cm~'^ for an un- 
shocked cloud with B = 5 ^iG and B = 10 ^.G of Fig. O 
respectively, before the interaction. 

As in Figure[T] the symbols (i.e., the crosses and the tri- 
angle) in Fig. [2] indicate the initial conditions assumed for 
the SNR-clouds interactions examined in the numerical sim- 
ulations described in Paper I for unmagnetized SNR-cloud 
interactions. We see that when the magnetic field is included, 
both the crosses and the triangle lie outside the SF domain 
of the diagrams. This means that for the initial conditions 
corresponding to these points in the diagram, SF is unlikely 
to occur. In paper I, the appfication of the results of the 
diagram of Figure [1] for an unmagnetized cloud with rc = 
10 pc and Uc ^ 10 cm~^ to the young stellar association 
of /3— Pictoris (see the region near the cross in the third 
panel of Fig. [1} had led us to conclude that this stellar as- 
sociation could have originated from recent past interaction 
between a cloud and an SNR with a radius of approximately 
52 pc. However, the inclusion of an effective magnetic field 
in the cloud with an intensity of only 1 ^.G has put the same 
cross corresponding to the initial conditions for formation of 
this young stellar association outside of the SF zone (see the 
cross in the third panel of Fig. [2]3. 



4 SNR-CLOUD INTERACTION: MHD 
NUMERICAL SIMULATIONS OF 
SELF-GRAVITATING CLOUDS 

As remarked before, there has been several numeri- 
cal studies of the impact of shock fronts on interstellar 




2000 4000 6000 8000 10000 



140 
120 

_100[i 
& 80 ■ 
d 60 
40 
20, 





// 


n // 






i \ ' 

















100 200 300 

nc (cm-^) 



400 



500 



150 

5:100 

CC 

z 

50 



W ' 




Wi 

; A 




\ \ 

■ +V- 














* We must remember, 
paragraph. 



however, the comment of the previous 



Figure 2. Constraints on the SNR radius versus cloud density 
for 4 different cloud radius as in Figure [ij but in the presence 
of a magnetized cloud with a normal B to the shock front B= 1 
/iG. Dashed (yellow) line: upper limit for complete cloud destruc- 
tion after an encounter with an adiabatic SNR derived from Eq. 
II25I I; solid (magenta) line: upper limit for the shocked cloud to 
reach the Jeans mass derived from Eq. I|17| l: dotted (blue) lines: 
upper limits for the shock front to travel into the cloud before 
being decelerated to subsonic velocities derived from Eq. I|30|l for 
different values of the cooling function A{Tgf^{tst)) = 10~^^ erg 
cm"^ s"-"^ (lower curve), 5 X lO"^'^ erg cm^ (middle curve), and 
3x 10~^^ erg cm"^ s~^ (upper curve). The dashed-dotted line (red) 
corresponds to the Jeans mass constraint for null magnetic field 
(see text for details). The shaded area defines the region where 
star formation can be induced by a SNR-cloud interaction (be- 
tween the solid, dashed and dotted lines). The crosses and the 
triangle in the panels indicate the initial conditions assumed for 
the clouds in the numerical simulations described in Section 4.2 
of Paper I. The star corresponds to the initial conditions 
of the MHD simulation of Figure [6) 
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Figure 3. The same as in Figure [2] for a cloud with Tc = 10 pc 
and now considering B = 5 fiG (upper panel) and B = 10 ^tG 
(bottom panel). 



clouds (e.g., Sgro 1975; McKee & Cowie 1975; Woodward 
1976; Nittmann et al. 1982; Tenorio-Tagle & Rozyczka 1986; 
Hartquist et al. 1986; Bedogni & Woodward 1990; Mac Low 
et al. 1994; Klein, McKee & Colella 1994; Anderson et al. 
1994; Dai & Woodward 1995; Xu & Stone 1995; Jun, Jones 
& Norman 1996; Redman, Williams & Dyson 1998; Jun & 
Jones 1999; Lim & Raga 1999; de Gouveia Dal Pino 1999; 
Miniati, Jones & Ryu 1999; Poludnenko, Frank & Black- 
man 2002; Fragile et al. 2004; StefTen & Lopez 2004, Raga, 
de Gouveia Dal Pino et al. 2002; Fragile et al. 2005; Marcol- 
ini et. al 2005; Melioli, de Gouveia Dal Pino & Raga 2005; 
Nakamura et al. 2006) most of which were mainly concerned 
with the effects of these interactions upon the destruction of 
the cloud. In particular, the studies that incorporated the 
effects of the radiative cooling revealed the relevance of it in 
delaying the destruction of the cloud and the mixing of its 
materials with the ISM (e.g., Melioli, de Gouveia Dal Pino 
& Raga 2005). 

In paper I, in order to check the predictions of our 
semi-analytic diagrams built for interactions involving SNR 
shocks and non-magnetized clouds, we performed 3-D hydro- 
dynamical radiative cooling simulations following the initial 
steps of these interactions. Here, we repeat this analysis but 
also take into account the effects of the magnetic fields and 
the self-gravity, in order to follow the late evolution of the 
shocked material within the magnetized clouds and check 
whether it suffers gravitational collapse or not, in consis- 
tence with our diagrams. 

To this aim, we have employed the grid code Godunov- 
MHD presented in Kowal & Lazarian (2007), and tested in 
Falceta-Gongalves et al. (2008), which solves the gas dynam- 
ical equations in conservative form as follows: 



dp 
dt 
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dB 
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V X (v X B) = 0, 



-pV*, (34) 



(35) 



+ V • (pv) = 0, 



(33) 



with V ■ B = 0, where p, v and p are the plasma density, ve- 
locity and pressure, respectively, B is the magnetic field and 
V'^^ — AirGp. The equations are solved using a second-order 
Godunov scheme, with an HLLD Riemann solver to properly 
consider the MHD characteristic speeds. For the self-gravity 
term, we used the FACR (Fourier Analysis Cyclic Reduc- 
tion) Poisson solver at each time step. In Paper I, given the 
importance of the radiative cooling behind the shocks, we 
simulated explicitly its effects in the hydrodynamical sim- 
ulations then presented. Presently, since the main focus in 
the simulations is to study the role of the magnetic field 
and self- gravity on the evolution of the clouds, we solved 
the equations under the approximation of strong radiative 
cooling. The set of equations above is closed by the equation 
of state p oc p"^, setting an effective 7e// = 1-2 to simulate 
the strong radiative cooling. 

The computational box has dimensions 100 pc x 100 pc 
X 100 pc, corresponding to a fixed mesh of 256'^ grid points. 
A SNR is generated by the explosion of a SN with energy 
Eq = 10^^ erg initially injected at the left-bottom corner 
of the box. Several runs were carried out with different ini- 
tially uniform cloud densities (ric), radii (rc), and distance 
between the SNR and the cloud's surface (Rsnr)- The ini- 
tial conditions are described in Table 1. We have selected 
two of these simulations to show in detail, as follows. 

Fig.|3]depicts the maps of the density and magnetic field 
vectors at the mid x-y plane that intercepts both the cloud 
and de SNR for a cloud with ric — 10 cm""^, Vc — 10 pc, 
Rsnr = 50 pc. These initial conditions represent the cross 
at the third panel of Fig. [2] (from top to bottom) which lies 
well outside of the SF shaded zone in the diagram and pre- 
dicts cloud destruction due to the impact. In the simulation 
of Fig. 4, the SNR shock front compresses the cloud increas- 
ing the magnetic energy density. After the interaction (at 
t ~ 1.8 X 10® s, a tail of gas is swept behind the cloud by 
the expanding SNR, similarly to the results obtained, 
e.g., by Murray et al. (1993). The magnetic energy den- 
sity along the tail is also increased. We see that, due to 
the compression, the cloud initially becomes gravitation- 
ally unstable and starts to collapse. At 1.8 x lO^yr, the core 
is ~ 10'' times denser. The contraction of the cloud also 
causes the increase of the magnetic energy that acts 
against further collapse. Later, the simulated cloud 
rebounds and undergoes expansion and evaporation, 
as seen at 3.3 x lO'^yr. This re-expansion may be oc- 
curring mostly because of the inefficient cooling of 
the simulation. Since we set 7e// = 1.2 the gas may be 
not being realistically cooled, as one should expect 
for a real ISM collapsing cloud. In such situations, 
7e// may be even smaller than unity, i.e. the shocked 
material may present temperatures lower than right 
before the shock. On the other hand, Spaans & Silk 
(2000) modeled the chemistry, thermal balance and 
radiative transfer for different conditions of the ISM, 
and showed that a polytropic pressure equation may 
be used as first approximation. The effective poly- 
tropic index is 7e// ~ 0.8—1.4, depending on the local 
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Table 1. Description of the simulations 



ric (cm (pc) Rsnr{p<^) code B {fJ.G) Result Prediction 



10 20 40 Hydro - collapse collapse 

10 20 40 MHD 1 stable stable 

50 5 15 MHD 1 stable stable 

10 10 50 MHD 1 stable stable 

100 10 25 MHD 1 collapse collapse 




X-AXIS X-AXIS 



Figure 4. Colour-scale maps of the mid-plane density distribution and magnetic field vectors evolution for the model with Uc = 10 
cm~'^, rc = 10 pc, RsNR = 50 pc, and B = 1/iG. The evolutionary times are t = 3 X lO^yr (top left-hand panel); 8 X lO^yr (top 
right-hand panel); 1.8 X lO^yr (bottom left-hand panel); and 3.3 X lO^yr (bottom right-hand panel). The ISM where the SNR expands 
has a number density n = 0.05 cm~^, and a temperature 10^ K. The field vectors are normalized by the maximum at each frame. 



conditions. Nonetheless, even though the real cloud 
would not re-expand in the presence of a stronger 
radiative cooling, the non-collapsing condition is still 
possible due to the increase in the magnetic energy. 
In any case, we are currently implementing a more 
realistic method to calculate the cooling function in 
our Godunov-MHD code, based on the interpolation 
method for a table of A(r) parameter (see Stone et 
al. (2008)). We intend to compare the presented re- 
sults with more realistic calculations in a forthcom- 
ing work. 

In Paper I, where non-equilibrium radiative 
cooling has been properly taken into account, but no 
magnetic field or self-gravity were considered, the 
hydrodynamical simulations (Fig. 4 of that paper) 
suggest that the cloud evaporates. As the magnetic 
field is introduced (Fig. [4] of the present paper) the 



increasing magnetic pressure at the later stages of 
the cloud evolution prevents the collapse, as pre- 
dicted by the SF diagram of Fig. [2j 

Figure [5] shows the distribution of several phys- 
ical parameters of the model of Fig. 4 at t = 1.8 x lO" 
yr. The column density is A'^ > 10^" cm~^ in the 
core. As remarked above, we see that the dense 
core formed during the initial collapse results an 
increase of the total gas and magnetic pressures 
with a maximum temperature T > 10"* K, while the 
magnetic energy density is an order of magnitude 
larger than that of the surrounding medium. The 
large increase of the gas pressure is probably the 
main responsible for the re-expansion of the cloud. 
This is highly dependent on the radiative cooling. 
In a more realistic situation, the cooling timescale 
tcooi ~ kT/nA(T) should be generally shorter than the 
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Figure 5. Colour-scale maps of the mid-plane density with magnetic field vectors, column density, gas pressure and magnetic energy 
density for the same model as in Figure |4] at t = 1.8 X lO^yr. 



collapsing timescale tg ~ {Gp)~^^, which means that 
the core temperature should be smaller than lO^K. 
Nevertheless, the magnetic field also plays a crucial 
role in stabilizing the cloud and supporting the cloud 
against collapse, as discussed below in Fig. [7J 

Figure |6] shows the mid-plane maps of the density distri- 
bution and the magnetic field vectors for the model of Table 
[T]with a cloud with ric = 100 cm~^ and Tc = 10 pc, and a 
SNR with RsNR ~ 25pc. The evolution in this case is similar 
to the previous model. However, after t — 1.2 x 10® yr, as the 
cloud collapses the magnetic pressure is not able to counter- 
balance gravity and the cloud keeps contracting. The higher 
cloud density in this case causes an efficient radiative cool- 
ing of the shock compressed material in the cloud that keeps 
the thermal energy low. Thus the collapse simply drags the 
magnetic field lines that increase the magnetic energy den- 
sity, but this will never be able to stop the collapse. This 
result is consistent with the predictions of the SF formation 
diagram. The initial conditions of this system of Fig. |6] cor- 
respond to the star symbol at the third panel of Fig. [2] (from 
top to bottom) which lies inside the SF shaded zone. 

As remarked before, the stability of a cloud supported 
by the magnetic pressure may be quantified by the mass-to- 
magnetic flux ratio, M/<1? ~ N/B, where N is the cloud col- 
umn density (e.g., Crutcher 1999). We have computed this 
mass-to-fiux ratio for the simulated clouds above at several 
time-steps. The results are shown in Figure 7. The open 
dots represent the maximum column density for each snap- 
shot as a function of the magnetic field averaged along the 



given line of sight (with maximum TV). We notice that for the 
model with Uc = 10 cm~^, = 10 pc and Rsnr ~ 50 pc (of 
Fig. 4), the cloud starts collapsing right after the interaction 
with the SNR. The column density increases and reaches the 
unstable regime. The cloud contracts and the total energy 
increases. Both the internal and the magnetic energy densi- 
ties suppress further collapse and then the cloud re-expands. 
Despite of the probably unrealistic re-expansion, as 
discussed before, the final stage is a stable cloud. For 
the case of ric — 100 cm^'', Tc ~ 10 pc, and Rsnr ~ 25 
pc, the same process occurs initially. However, the 
increase of the magnetic and internal energy densi- 
ties are not enough to avoid the continuous collapse. 

The other simulated models of Table [1] have presented 
results which are also consistent with the SF diagrams. The 
first model of the table is a pure hydrodynamical system 
with no magnetic field whose initial conditions correspond 
to the cross labeled in the bottom panel of Fig. [1] According 
to Table [TJ the numerical simulation of the evolution of this 
SNR-cloud interaction including the effects of self-gravity 
leads to the gravitational collapse of the cloud, which is in 
agreement with the prediction of the SF diagram. When a 
magnetic field of 1 /iG is included in a system with the same 
initial conditions, the MHD numerical simulation shows that 
the magnetic pressure prevents the cloud collapse (second 
model of Table [l| . This model is represented by the cross 
in the bottom panel of Fig. [2] which consistently lies outside 
the gravitational unstable zone of the diagram. 

The third model of Table[T]also lies outside the unstable 
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Figure 6. Colour-scale maps of the mid-plane density distribution and magnetic field vectors for the model with ric = 100 cm~ , = 10 
pc, RsNR = 25 pc, and B = 1/iG. The evolutionary times are t = 2 X 10^ yr (top left-hand frame); 5 X 10^ yr (top right-hand frame); 
1.2 X lO^yr (bottom left-hand frame); and 2.5 X lO" yr (bottom right-hand frame). The vectors are normalized by the maximum at each 
frame. 
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Figure 7. Mass-to-flux relation obtained for several time-steps of the simulations described in the text, for cloud conditions n = 10 
cm~^, Tc = 10 pc, RsN = 50 pc (left), and n = 100 cm"-^, Tc = 10 pc, Rsn = 25 pc (right). The solid line in each diagram gives 
the critical mass to flux ratio which separates the gravitational unstable (on the right-hand side of the line) from the stable domain 
(loft-hand side of the line). 



shaded zone of the SF diagram (see the X symbol in the 
second panel of Fig. and the MHD numerical simulation 
of the evolution of this system shows that the cloud ends up 
evaporating due to the strong SNR shock interaction. 



5 APPLICATION TO ISOLATED REGIONS OF 
THE ISM 

We can apply the simple analytical study above to isolated 
star formation regions of our own galaxy. In paper I, we 
focused on the formation of the young stellar association of 
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/3— Pictoris induced by a SNR-cloud interaction. Here, we 
will address few other examples in our ISM that present 
some evidence of recent past interactions with SNRs, like 
the Large CO Shell in the direction of Cassiopea (Reynoso 
& Mangum 2001) and the so called Edge Cloud 2 in the 
direction of Scorpions (Ruffle et al. 2007, Yasui et al. 2006) . 
A counter example is the region apparently without star 
formation around the Vela SNR. In fact, we will see below 
that the conditions of this region correspond to a point in 
our diagrams that lies outside the shaded SF zone. 

The Large CO shell is an expanding structure with a 
velocity ~ 3 km/s, a mass of 9.3 x 10^ Mq and a den- 
sity of ~ 35 cm~^. Reynoso & Mangum (2001) suggest that 
this expanding structure has probably originated from the 
explosion of a SN about ~ 4 x 10^ yr ago. Assuming that 
the cloud mass was originally uniformly distributed within a 
sphere of radius of 435" (~ 50 pc), the initial density could 
be Uc — 30 cm~'^. The SN shock front possibly induced the 
formation of the O 9.5 type star that has been detected as 
an IR source (IRAS 17146-3723). The Large CO Shell has 
an external radius of 50 pc and an inner radius of ~ 28 pc 
(Reynoso & Mangum, 2001). The age and small expansion 
velocity suggest that it is now a fainting evolved SNR. If 
we consider a cloud with the above density and radius (~ 
50 pc) at the time of the potential interaction with SNR in 
the adiabatic regime, we can identify this system in the SF 
diagram within the shaded zone, as indicated in Figure [S] 
if the SNR had a radius between 2.5 — 72 pc and an am- 
bient medium density n ~ 1 cm^'^. However, when we 
include a magnetic field in the cloud of 1 fiG, the range of 
possible radii for the SNR is reduced to Rsnr ~ 7.8 — 72 pc 
if the maximum radius is calculated using a radiative cool- 
ing function A = 3 x lO"^'^ erg cm"^ s~^. For an average A 
= 5 X 10~^^ erg cm'' s~^, the maximum possible radius of the 
SNR is reduced to 38.4 pc. Considering that the present ra- 
dius of the evolved SNR is probably around 50 pc, the range 
above of initial conditions for the SNR-cloud interaction is 
quite plausible. 

The SNR Vela has an almost spherical, thin HI shell 
expanding at a velocity of ~ 30 km/s. Instead of imping- 
ing on an interstellar cloud, it is expanding in a fairly dense 
environment with evidence of some structure formation. As- 
suming that Vela is at a distance of ~ 350 pc from the Sun, 
its shell radius is of the order of 22 pc. The ambient density 
is ~ 1 to 2 cm~^ and the initial energy of the SN was around 

1 - 2.5 X 10^^ erg (Dubner et al. 1998). These initial con- 
ditions correspond to the square symbol in the diagram of 
Figure |8] and it lies outside the SF shaded zone. This is con- 
sistent with the absence of dense clouds, clumps, filaments, 
or new born stars in the neighborhood of this SNR. 

The SF region in the neighborhood of the Edge Cloud 

2 is a possible example involving an interaction of a SNR 
in the radiative phase with a cloud. This is actually a giant 
molecular cloud with a diameter ~ 30 — 40 pc. It is one of 
the most distant cloud complexes from the galactic center of 
the Milky Way (d ~ 22 - 28 kpc; Ruffle et al. 2007) . For this 
reason, it is embedded in a region where the gas pressure 
is extremely small and the presence of external SF agents 
like spiral arm perturbations is improbable. This could be 
an indication that this cloud complex is very stable, except 
for the recent detection of two associations of T-Tauri stars 
with ages of ~ 10^ yr (Kobayashi & Tokunaga 2000, Yasui 
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Figure 8. Diagram presenting the parameter space domain for 
SF for a cloud interaction with a SNR in the adiabatic regime. 
The light-gray shaded zone corresponds to the allowed SF zone 
for the interaction without magnetic field while the dark-gray 
corresponds to the interaction with the magnetized cloud with 
B = 1 fiG. As in the previous figures, the dotted (blue) lines 
from top to bottom represent the shock penetration constraint 
for three different values of A = 3 X 10~^^ erg cm^ s~^; 5 X 
10~^^ erg cm^ s~^; and 1 X 10~^^ erg cm-^ s~^, respectively. 
The other initial conditions are rc = 50 pc and Tc = 100 K. 
The ambient medium density is n = 1 cm~'^. The triangle 
represents the average conditions for the interaction involving the 
Large CO Shell, centered at ric = 30 cm~^ and Rsnr = 36 
pc, while the (red) ellipse indicates a possible range of values 
for the interacting SNR radius. The square represents the initial 
conditions for Vela system (see text for details). 



et al. 2006). The Edge 2 cloud has a temperature of 20 K, 
a density n/^j ~ 10* cm~^, and an estimated mass of ~ 
10" Mq. There is an old and large SNR associated to this 
cloud, GSH 138-01-94, which consists of an HI shell with a 
radius of 180 pc expanding into the ISM with a speed of 
11.8 ± 0.9 km/s and an age of about 4.3 x 10® yr (Ruffle et 
al. 2007). According to Ruffle et al. (2007), the formation of 
the present structure and chemical composition of the Edge 
2 Cloud is possibly a result of interactions of this SNR with 
the IS gas. 

Considering the SNR characteristics above (i.e., its ve- 
locity, radius and an energy ~ 10^^ erg) and using Eq. 8 of 
the Paper I, we can estimate an ambient density n ~ 0.14 
cm~'^. Also, assuming that the mass of the cloud was orig- 
inally distributed uniformly in a sphere of average radius 
^ 17.5 pc, we find that the density of the cloud before the 
compression was Uc ~ 14 cm~^. With these initial condi- 
tions, the SNR-cloud interaction would lie within the al- 
lowed SF zone of the diagram of Figure [9] for a SNR with a 
radius between ~ 31 and 102 pc at the time of the interac- 
tion and for a value of A = 3 x 10~^^erg cm'' s~^ . We can try 
to constrain further the possible value of the radius of the 
SNR assuming that the interaction started sometime after 
the remnant had become radiative and before 10® yr, which 
is the approximate age of the observed stars. These limits 
imply Rsnr — 46 pc - 84 pc, as indicated by the region 
within the ellipse in the diagram of Figure (9) 
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Figure 9. Diagram presenting the parameter space domain for 
SF for a cloud interaction with a SNR in the radiative regime 
expanding in an ambient medium with n ~ 0.14 cm~^. The light- 
gray shaded zone corresponds to the allowed SF zone for the in- 
teraction without magnetic field while the dark-gray corresponds 
to the allowed SF zone for the interaction with the magnetized 
cloud with B = 1 fiG. As in the previous Figures, the dotted 
(blue) lines from top to bottom represent the shock penetration 
constraint for three different values of A = 3 X 10 "■^^ erg cm'^ s~^ ; 
5 X lO"'^^ erg cm^ s~^; and 1 X 10~^^ erg cm^ s~^, respectively. 
The other initial conditions are rc = 20 pc and Tc = 100 K which 
are appropriate for the Cloud 2-SNR GSH 138-01-94 system. 
The ellipse indicates a possible range in the parameter space for 
the interaction between the Edge Cloud 2 and the SNR to result 
in gravitational collapse. It spans the ranges ric ~ 9.4 — 22 cm^'^ 
and RsNR ~ 46 — 84 pc. 



6 ESTIMATING THE EFFICIENCY OF STAR 
FORMATION FROM SNR-CLOUD 
INTERACTIONS 

In the study we have carried out here, we concentrated on 
isolated interactions between diffuse clouds and SNRs with- 
out focusing on the effects that such interactions can have 
upon the global SF in the Galaxy. The present star forma- 
tion efficiency is typically observed to be very small, of the 
order of a few 0.01 in dispersed regions, but it can attain a 
maximum of ~ 0.3 in cluster-forming regions (Lada & Lada 
2003; see also Nakamura & Li 2006 for a review). 

From the present analysis, we can try to estimate the 
star formation efficiency that interactions between SNRs 
and diffuse clouds produce and compare with the observed 
values in order to see the contribution of this mechanism on 
the overall SFE in the Galaxy. The diagrams built in this 
work provide a SF domain for these interactions, in other 
words, they establish a set of conditions that isolated inter- 
actions must satisfy in order to result in successful gravita- 
tional collapse of the compressed cloud material. In order to 
evaluate the corresponding global SFE of these interactions 
we have to calculate first their probability of occurrence in 
the Galaxy. Considering that once a SNR is formed it will 
propagate and compress the diffuse medium around it pro- 
viding the sort of interactions we are examining, the prob- 
ability of these interactions to occur must be proportional 
to: 



NsNII'^tsNRiRs 



. Asnr{Rs 



10" 



Ats 



1.7 X 10-2s 5 X 105s 2.5 x lO'^pc^ 4 x 10»pc' 



(36) 
(37) 



Where we have assumed a homogeneous galactic thin disk 
with a radius of 20 kpc to compute the galactic area {Aq), 
and where Nsnii is the rate of SNII explosions (e.g. Cap- 
pellaro, Evans & Turatto 1999), AtsNniRsNR) is the life- 
time and AsNR is the area of a SNR and both depend on 
RsNR- Since not all the galactic volume is filled up with 
clouds, this quantity above must be multiplied by the dif- 
fuse neutral clouds filling factor in order to give the ap- 
proximate probability of occurrence of SNR-clouds interac- 
tions. If we consider the amount of gas that is concentrated 
within cloud complexes in the cold phase of the ISM, the 
corresponding volume filling factor of the clouds is /c — 5% 
(e.g., de Avillez & Breitschwerdt 2005), and then the prob- 
ability of occurrence of the interactions will be given by 
fsNR-c fsNR X /c. Substituting Eq. (2) and (6) for 
A.tsNR{RsNR) of Paper I in Eq. P?)) we find: 

fsNR-c{RsNR,a) = fsNR ^ fc — 9.5 X 10~^ Rs^NR,50 

for an interaction with a SNR in the adiabatic regime and 

fsNR-c{RsNR,r) = 2.9 X 10"^ RfNR,50 

for an interaction with a SNR in the radiative regime. 

This probability can be multiplied by the mass fraction 
of the shocked gas that is gravitationally unstable within 
the SF domain of our SNR-cloud interaction diagrams in 
order to give an effective global star formation efficiency for 
these interactions. Using the calculated Jeans mass for the 
shocked gas (Eqs. 1151 \W\ and I20|l as an approximate lower 
limit for the mass fraction of the cloud that should collapse 
to form stars, we obtain: 



sfesNR-c — fsNR- 



mj 



which in the case of a non-magnetized cloud interacting with 
a SNR in the adiabatic regime gives (according to Eq. I15|) : 

rn2 Jj6 



sfesp 



5.4 X 10" 



Tc,l00 -n-SJVfl.SO 

nc,io ri^o E51 h 



(38) 



And for the case of the interaction of a non-magnetized 
cloud with a SNR in the radiative phase (using Eq. I16p : 



sfesNR 



1.5 X 10" 



-'c.lOO ^SNR.50 " 



(39) 



For the interaction involving a magnetized cloud we 



have; 



afesNR-c,B 



1/2 „3 
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where y must be substituted by Eq. [S] Using Eqs. |9] and 
as] it gives the sfe for interactions with adiabatic SNRs and 
using Eqs. [10] and IA9I it gives the sfe for interactions with 
radiative SNRs. 

As examples. Figures [10] and [11] show plots of the ap- 
proximate sfe computed for SNR-clouds interactions as a 
function of the SNR radius for different values of the cloud 
magnetic field and for different values of the cloud density. 
Since the computed sfe is already dependent of the Jeans 



mass constraint (see Egs. 1381 and I39|l . the zone that defines 
the allowed SF domain in these new plots is constrained by 
the other two conditions over Rsnr, i.e., the shock penetra- 
tion into the cloud and the cloud non-destruction conditions 
calculated in Section 3. The increase of the magnetic field 
tends to shift the SF zone to smaller values of the radius of 
the SNR both for SNRs in the adiabatic and in the radiative 
regimeslfl The vertical line in the diagrams establishes the 
transition radius of the SNR from the adiabatic (left) to the 
radiative (right) regime for the specific initial conditions of 
the diagrams. In other words, the relevant domain in the 
diagrams of Figure [10] is on the left-hand side of this line, 
while in Fig. 1111 it is on the right-hand side. 

We note that the effective sfe found for these interac- 
tions is generally smaller than the typical values observed 
for the Galaxy and the range of possible values for sfe de- 
creases with the increase of the magnetic field in the cloud. 
We also see that the parameter space that allows sfe values 
close to the observed values (0.01—0.3) is very reduced. The 
diagrams of Fig. 1111 which appy to interactions involving 
more evolved SNRs already in the radiative phase, shows 
interesting results. In this phase the SNR is less powerful 
and therefore, less destructive than in the adiabatic phase. 
Besides, it is much more expanded increasing the probability 
of encounters with clouds. This explains the larger attained 
values of the calculated sfesNR-c- However, the parame- 
ter space that allows successful interactions for SF on the 
right-hand side of the diagrams is even thinner than in the 
interactions with adiabatic SNRs (Fig. llOf) and disappears 
when Be = 10 fj.G. 

The results above suggest that these interactions are 
not sufficient to explain the observed sfe of the Galaxy ei- 
ther in the presence or in the absence of the magnetic field in 
the cloud. They are consistent with previous analysis per- 
formed by Joung & Mac Low (2006) where these authors 
have concluded that Supernova-driven turbulence tends to 
inhibit global star formation rather than triggering it. We 
should note however, that they have based their conclusion 
on the computation of the star formation rate (SFR), rather 
than the sfe, from box simulations of the ISM with SN tur- 
bulence injection and their computed SFR has been weighed 
by a fixed value of the sfe taken from the observations (sfe 
~ 0.3). 



7 CONCLUSIONS 

We have presented here a study of isolated interactions be- 
tween SNRs and difi'use neutral clouds focusing on the deter- 
mination of the conditions that these interactions must sat- 
isfy in order to lead to gravitational collapse of the shocked 
cloud material and to star formation, rather than to cloud 



^ As remarked in Sec. 3.4, this shift that is imprinted 
in the Jeans mass constraint when considering a non- 
null magnetic field normal to the shock front is expected 
to affect only the initial compression and collapse of the 
shocked cloud since the later collapse will occur mostly in 
the direction parallel to B. Thus, in Figs. 10 and 11, the 
top diagrams should be considered as the most realistic 
ones when comparing with the observations. 
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Figure 10. Plots of the calculated star formation efficiency as 
a function of the SNR radius in the adiabatic regime for several 
values of the cloud density represented by the continuous (red- 
dish) lines (from top to bottom: 10, 20, 40, 60, 80, 100, 150, 200, 
300, 400, 500, 600, 700, 800, 900 cm-^') for an interaction with 
a cloud with: (a) Be = (top panel); (b) Be = 1 /iG ( middle 
panel); and (c) Be = 10 /zG (bottom panel). The other initial 
conditions are = 10 pc, Te = 100 K, and A = 5 X 10~^^ erg 
cm"^ s"'^. The ambient medium density used is n = 0.05 
cm"''. The two constraints, i.e., the non-destruction condition 
of the cloud (solid, green line) and the shock penetration condi- 
tion (dotted, blue line) delineate the zone for allowed SF in these 
sfe versus Rsnr diagrams. The continuous vertical (black) line 
represents the transition radius of the SNR from the adiabatic 
to the radiative regime so that only the values of Rsnr on the 
left-hand side of this line are relevant in these diagrams (see text 
for details). 

destruction. A preliminary study of these interactions ne- 
glecting the effects of the magnetic field in the cloud had 
been performed previously by Melioli et al. (2006) . We have 
presently incorporated these effects and derived a new set of 
conditions for the interactions. A first condition determines 
the Jeans mass limit for the compressed cloud material in 
the impact. A second one establishes the penetration ex- 
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Figure 11. The same as in Figure [TOl but considering the inter- 
action with a SNR in the radiative phase. In these diagrams, only 
the values of Rsnr on the right-hand side of the vertical (black) 
line are relevant in this phase. 



tent of the SNR shock front inside the cloud before being 
stalled due to radiative losses. It must have energy enough to 
compress as much cloud material as possible before fainting. 
A third constraint establishes the condition upon the same 
shock front under which it will not become too strong to 
destroy the cloud completely. We have then built diagrams 
of the radius of the SNR as a function of the cloud den- 
sity where this set of constraints delineate a domain within 
which star formation may result from these SNR-cloud in- 
teractions (Section 3). As expected, we find that an embed- 
ded magnetic field in the cloud normal to the shock front 
with an intensity of 1 fiG inhibits slightly the domain of SF 
in the diagram when compared to the non-magnetized case. 
The magnetic field plays a dominant role over the Jeans con- 
straint causing a drift of the allowed SF zone to higher cloud 
densities in the diagram. When larger intensities of magnetic 
fields are considered (5-10 fJ,G), the shrinking of the allowed 
SF zone in the diagrams is much more significant. We must 
emphasize however that, though observations indicate typ- 



ical values of Be — 5 — lOpG for these neutral clouds, the 
fact that we have assumed uniform, normal fields in the in- 
teractions have maximized their effects against gravitational 
collapse. We should thus consider as more realistic the result 
obtained when an effective Be — IfiG was employed. These 
diagrams derived from simple analytical considerations pro- 
vide a useful tool for identifying sites where star formation 
could be triggered by the impact of a SN blast wave. 

We have also performed fully 3D MHD numeri- 
cal simulations of the impact between a SNR and a 
self-gravitating cloud for different initial conditions 
(Section 4) tracking the evolution of these interac- 
tions and identifying the conditions that have led ei- 
ther to cloud collapse and star formation or to com- 
plete cloud destruction and mixing with the ambient 
medium. We have found the numerical results to be 
consistent with those established by the SNR-cloud 
density diagrams, in spite of the fact that the later 
have been derived from simplified analytic theory. 
We remark that the radiative cooling in the MHD 
numerical simulations has been considered through 
the adoption of an approximate polytropic pressure 
equation with 7e// ~ 1.2 (see Spaans & Silk 2000). 
We are presently implementing a more realistic cool- 
ing function in our Godunov-MHD code in order to 
test, e.g., the re-expanding cloud case and its valid- 
ity under more realistic situations. 

We have applied the results above to a few examples of 
regions in the ISM with some evidence of interactions of the 
sort examined in this work. In paper I, the application of the 
results of the SF diagram for a non-magnetized cloud (as in 
Figure [2} to the conditions around the young stellar associ- 
ation of /3— Pictoris in our local ISM had led us to conclude 
that this stellar association could have originated from past 
cloud-SNR interaction only under very restrict conditions, 

1. e., with a cloud with radius ~ 10 pc and density ~ 20 cm""^ 
and a SNR with a radius ~ 42 pc. However, in the present 
work we find that with the inclusion of an effective magnetic 
field in the cloud with an intensity of only 1 /iG this inter- 
action is unlikely to produce that stellar association (Figure 

2, cross in the third panel from top to bottom), at least not 
for the set of initial conditions proposed in the literature for 
that system (see also Melioli et al. 2006). In the case of the 
expanding Great CO Shell— 09. 5 star system, we find that 
local star formation could have been induced in this region 
if, at the time of the interaction, the SNR that probably orig- 
inated this expanding shell was still in the adiabatic phase 
and a radius between ~ 8 pc — 29 pc impinged a magnetized 
cloud with density around 30 cm"'' (Figure (Sjl. Another ex- 
ample is the SF region near the Edge Cloud 2. This is one of 
the most distant cloud complexes from the galactic center 
where external perturbations should thus be rare. But the 
recent detection of two young associations of T-Tauri stars 
in this region could have been formed from the interaction of 
a SNR in the radiative phase with a cloud, if the interaction 
started < lO" yr, the SNR had a radius Rsnr ~ 46 pc — 84 
pc and the magnetized cloud a density around Uc ~ 14 cm""' 
(Figure [9ll. 

Finally, though in this study we have focused on isolated 
interactions involving SNRs and clouds, we used the results 
of the diagrams to estimate the contribution of these inter- 
actions to global star formation. Our evaluated effective star 
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formation efficiency for tliis sort of interaction is generally 
smaller than the observed values in our own Galaxy (sfe ~ 
O.Of — 0.3) (Figures [To] and [TT]). This result seems to be con- 
sistent with previous analysis (e.g., Joung & Mac Low 2006) 
and suggests that these interactions are powerful enough to 
drive structure formation, supersonic turbulence (see, e.g., 
simulation of Figure S} and eventually " local" star forma- 
tion, but they do not seem to be sufficient to drive global 
star formation in our galaxy or in other normal star forming 
galaxies, not even when the magnetic field in the neutral 
cloud is neglected. In conclusion, the small size of the al- 
lowed SF domain in the diagrams and the results for the es- 
timated sfe indicate that these interactions must lead more 
frequently to the destruction of the clouds, rather than to 
their gravitational collapse. 
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Figure Al. Schematic diagram of the interaction between an 
SNR and a cloud. The SNR expands and impacts the cloud. The 
angles /3 and 7 are functions of the time, the SNR velocity, and 
the cloud and SNR radii, as indicated by the equations of the 
text. 



APPENDIX A: 

For spherical SNR— cloud interactions we need to consider 
the effects of curvature in the cloud-SNR interactions. The 
instantaneous velocity of the shocked gas moving towards 
the center of the cloud, Vcs, is only a fraction of the SNR 
velocity and depends on the density contrast, Xi between 
the shell and the cloud (as in the planar shock case) and the 
angle 7 between the SNR velocity vector and the line that 
links the center of the cloud and the instantaneous contact 
point between the cloud and the SNR (Fig. IA1|I . We see 
that at t = 0, i.e, when the SNR touches the cloud, these 
two lines are coincident and 7 = 0, then Vcs = X'''^'fSNR- 
Later, at a time tc,sNR when the SNR approaches 
the center of the cloud and the SNR shock energy 
input ends, 7 = 7r/2 and Vcs = The average value of 
the velocity integrated over this SNR crossing time, is 



VSNR 



\ He J 



tc,SNR 



where 



cos7(t) 



RsNR{t) 



cos7(t) dt (Af) 



(A2) 



2RsNR{t) 

d is the distance between the SNR center and the cloud 
center, Rsnr is the SNR radius, is the cloud radius, ric 
is the initial cloud density, n^h is the SNR shell density and 
VSNR is the expansion velocity of the SNR. 

^ We note that only for a planar shock, 7 will be equal 
to 7r/2 when the shock reaches exactly the center of the 



cloud. For rc ^ RsNR, 7 will be 
reaches the center. 



tt/2 a little before it 



We assume that the distance between the center 
of the cloud and that of SNR remains constant: 

d = Rsnr + (A3) 

After a time tc,sNR, this distance can be written as 

d = y^RsNR{t)^ + r'i (A4) 

At this time, the SNR radius is approximately: 

RsNR{t) — Rsnr + vsnr tc,sNR ■ 
Thufl: 

V Rsnr + '^^'^ ^snr ~ Rsnr 

tc.SNR ^ , (A5) 

VSNR 

where 

Rsnr ~ Rsnr{0) . 

With the integration limits above, eq. (|A1|I has a solu- 
tion that should replace the more approximate one given 
in Paper 10 

Figure IA2I depicts the plot of the integral term 
of eq. (IATI) . I = -r^ — f*c.s]VH ^j. ^.j^^^j 

^ ^ " tc.SNR JO /V / I 

with rc — 10 pc as a function of Rsnr- It clearly 
shows that the effect of curvature will be relevant 
only for values of Rsnr/tc near unity. The corrections 
in the solution above to Vcs will produce a few changes in 

^ We note that in Paper I, it was assumed that the shock 
always reaches the center so that t^s SNR = ^c/vsnRj The 
factor 2 that appears below eq. (12) in Paper I is a typo. 

We note that the multiplying factor 2 to the integrand 
of eq. (lAlll . which is due to the fact that the total aper- 
ture angle of the SNR shock is actually 2/3, rather than 
(3, was not considered in Paper I. 
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Figure A2. Values of I for a cloud with rc = 10 pc and for 
different values of RsNR 



the multiplying factors that appear in Equations (13) to (26) 
of Paper I (in the absence of magnetic field), as indicated 
below: 
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where tec is the cloud crushing time 
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i.e., the time the internal shock takes to cross the 
cloud, M = Vca/'y^^^Ca is the Mach number of the shock 
into the cloud and Cs is the cloud sound speed, ric^sh is the 
shocked cloud gas density, rc,io is the cloud radius in units 
of 10 pc. Tic, 10 is the unshocked cloud density in units of 10 
cm""', £51 is the SN energy in units of 10^^ erg, Rsnr,50 
is the SNR shell radius in units of 50 pc, fio is the den- 
sity contrast between the SNR shell and the ISM density in 
units of 10, n is the ambient medium density, Tc,ioo is the 
cloud temperature in units of 100 K, and Is is the I factor 
calculated for Rsnr/tc = 5. In the equations above the in- 
dices "a" and "r" refer to interactions involving SNRs in the 
adiabatic and in the radiative phase, respectively. 
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